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1 Introduction 

During the last decade, the investigation of quark matter at finite temper- 
ature and density has become one of the QCD focal points. It is expected 
• r- 1 ! that at densities which are 3-5 times larger than the normal nuclear density 

baryons are crushed into quarks. It is also expected that if the temperature 
of such quark matter is low enough (below few tens of MeV) the system is un- 
stable with respect to the formation of quark-quark Cooper-pair condensate 
PQIH] This phenomenon is called color superconductivity since diquarks 
belong to the 3 color channel. At present there is a fair understanding of 
color superconductivity physics in the regime of ultra-high density when a s 
is small [4J. Most interesting is, however, the region of moderate densities 
(3-5 times larger than the normal nuclear density). This region is important 
for physics of neutron stars and may possibly be investigated in the labo- 
ratory, in particular, in future experiments at the GSI heavy-ion machine 1 . 

In heavy-ion collisions high-density state is naturally formed with high temperature 
thus preventing diquark condensation. 
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In the moderate density /strong coupling regime the theory faces the well- 
known difficulties of the nonperturbative QCD. Lattice QCD calculations 
encounter serious obstacle at nonzero density since in this case the determi- 
nant of the Dirac operator is complex resulting in non-positive measure of 
the corresponding path integral. Still, several attempts to perform lattice 
calculations at nonzero density have been performed - see the review paper 
0. The region of moderate densities has been extensively studied within 
the framework of the Nambu-Jona-Lasinio (NJL) type models (see d and 
references therein), using the instanton gas model [8J, or chiral perturbation 
theory 0. 

The main conclusion reached within the NJL-type models is that transi- 
tion from the nuclear matter (NM) phase to the quark matter (QM) phase 
occurs very early, namely when the quark density reaches the value only three 
times larger than the density of quarks in normal nuclear matter. This cor- 
responds to the value \x ~ 400 MeV of the quark chemical potential. In the 
NJL-type models the (NM) — > (QM) transition has two main signatures, 
namely: 

(i) The gap equation for the diquark channel acquires a nontrivial solu- 
tion. 

(ii) The quark constituent mass tends to zero 2 and the chiral symmetry 
is restored. At T = the transition is believed to be of the first order. 

Both conclusions should be taken with reservations due to oversimplifica- 
tions inherent for the models and certain arbitrariness in the interpretation 
of the results. 

The principal deficiency of the existing approaches to the (NM) — ► (QM) 
transition is the lack of understanding what happens to the gluon sector when 
the density increases. From thermodynamic arguments it follows that pos- 
sible color diquark condensate and gluon condensate have the same energy 
scale which results in competition between them somewhat similar to Meiss- 
ner effect [THj- It also has been demonstrated within the Ginzburg-Landau 
approach (see below) that fluctuations of the gluon field are important and 
may significantly influence the character of the phase transition and the criti- 
cal temperature [TJ [TTJ [J2] . From the fact that the would-be diquark conden- 
sate belongs to color anti-triplet it follows that 5 of 8 gluons become massive 
[3]. However, first principle derivation of the QCD string and gluon conden- 

2 In the density region under consideration only u- and d-quarks participate in possible 
pairing 



2 



sate evolution with the increase of the density is not currently available. An 
investigation aimed at the resolution of these problems has been attempted 
very recently [T5] . 

In the present paper we specify the set of the key parameters which char- 
acterize the (NM) — > (QM) transition and estimate the values of these 
parameters. In this way we obtain a model-independent though schematic 
picture which exhibits several nontrivial phenomena. The onset of the quark 
phase and its further evolution to higher densities may be viewed as a 
crossover from the strong coupling regime of composite nonoverlapping bosons 
(diquarks) to the weak coupling regime of macroscopic overlapping Cooper 
pair condensate 3 . Another feature of the transition region is the drastic in- 
crease of fluctuations. The (NM) — > (QM) evolution possesses also some 
features of the Anderson transition. All these three properties are interre- 
lated. 



2 BEC-BCS Crossover and the 
Ginzburg-Levanyuk Number 

As already mentioned in the Introduction, the gap equation for the diquark 
channel derived within the framework of the NJL-type models acquires a 
nontrivial solution (a gap) at \i ~ 0.4 GeV jHllZj- To describe the crossover 
we will need the corresponding value of the quark number density. Relation 
between \x and n is given by the well-known equation 

where f2(T, /i, A) is the thermodynamic potential, A - the gap parameter. For 
the NJL-type of models fl is easily calculated Ej (see below). According 
to jB] for Nf = 2 and T = transition to A ^ phase occurs at \i — 0.292 
GeV 4 . Then Eq. (0) yields n 1/3 = 0.18 GeV This number was obtained 
in the chiral limit , i.e. under the condition that the quark constituent mass 
goes to zero when A^O. In our view the conclusion that m = when A / 

3 This is a particular case of the Bose-Einstein condensate to Bardeen-Cooper-Schrieffer 
(BEC-BCS) crossover - see below. 

4 Strange quark starts to participate in pairing an much higher densities when ^ 3> 
TO, ^ 150 MeV. 
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may be a specific feature of the NJL model (or some versions of this model). 
The behavior of the quark constituent mass throughout the (NM) — > (QM) 
transition is still a moot point. 

The above result for n may be with a good accuracy reproduced using 
the equation for free degenerate quarks 

n = N c N f ^- 2 = I4 ~ V> (2) 

where \x = (kp + m 2 ) , and fj, — kp in the chiral limit. For \x = 0.292 GeV 
Eq.(j2J) yields n 1 / 3 = 0.17 GeV in a good agreement with Eq.ljT]). For fi = 0.4 
GeV Eq.© gives n 1/3 = 0.23 GeV. Thus we conclude that n 1/3 ~ 0.2 GeV 
~ 1 fm -1 in the transition region. 

We now turn our attention to the physics behind the fact that the gap 
equation derived within the framework of the NJL model and in the mean- 
field approximation acquires a nontrivial solution starting from \x ~ 0.4 GeV. 
It took quite some time before it was realized [T4*l ITS'] that a nonzero value 
of the gap does not mean the onset of the color superconductivity (the BCS 
regime). It is only a signal of the presence of fermion pairs. Depending on 
the dynamics of the system, on the fermion density, and on the temperature, 
such pairs may be either stable, or fluctuating in time, may form a BCS 
condensate, or a dilute Bose gas, or undergo a Bose-Einstein (BE) conden- 
sation. The fact that there is a continuous transition (crossover) from the 
strong coupling/low density regime of independent bound state formation to 
the weak coupling/high density cooperative Cooper pairing is well known [IH]. 
In contrast to macroscopic Cooper pairs, the compact molecular - like states 
which are formed in strong coupling/low density regime are called Schafroth 
pairs [HI HE]. 

The dimensionless crossover parameter is n 1//3 £ , where £ is the charac- 
teristic length of pair correlation when the system is in the BCS regime and 
the root of the mean square radius of the bound state when the system is 
in the strong coupling regime. Some arbitrariness occurs in the definition of 
the crossover parameter. For example, in [T^j it is defined as fci?£ ~ 1.7n 1//3 £ 
(see (J2J)). Another definition of the crossover parameter is xq = ///A |2"U] . 

In the BCS theory £ is given by £ = vp/irA, where vp is the velocity at 
the Fermi surface. For a typical metal superconductor Vp — c/137, A ~ 5K, 
so that £ ~ 10~ 4 cm. The density of electrons is n ~ 10 22 cm -3 . Therefore 
in the BCS regime n 1 ^ 3 ^ > 10 3 . In coordinate space the wave function of 
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the Cooper pair is proportional to (sin kpr/kpr) exp(— r/£) and hence it has 
~ 10 3 nodes. 

The crossover from the BCS to the strong coupling regime occurs at ^H] > 



The width of the crossover region with respect to the above parameter is 
several units and is model-dependent [IH]-[22]- It was first pointed out in 
[HI H3] that at ii ~ (0.3 — 0.5) GeV the quark system is in the crossover 
regime and not in the BCS regime as it was inferred from the fact that at 
such values of \x the gap equation acquires a nontrivial solution. 

Let us estimate the value of the crossover parameter at the onset of the 
phase with A ^ 0, i.e., at fi ~ (0.3 — 0.5) GeV. We have seen that n 1 / 3 ~ 
1 fm _1 in this region 5 . The value of £ in the strong coupling regime cannot 
be evaluated from the first principles. One may expect that it is of a typical 
hadronic scale £ ~ (1 — 2) fm and that it grows with density asymptotically 
approaching the BCS value. Model calculations confirm this expectations 
[2Sj "EHl- At zero density the root-mean square radius of the di quark is 
~ 1 fm [22] • At low density the single-gluon exchange model leads to the 
result £ ~ n^ 1 / 3 [23], while at fj, ~ 1 GeV ^n 1 ^ ~ 10 [25]. It should be noted 
that the calculation of the pair size in the crossover region is a complicated 
task even for the "simpler" system of ultracold fermionic atoms [2*^]. To sum 
up, we conclude that ^n 1 / 3 ~ 1 at the onset of the A ^ phase. The 
(NM) —>■ (QM) transition brings the system to the quark matter state in 
the crossover regime, but not in color superconducting state. 

Next we turn to the structure of the diquark wave function in the strong 
coupling regime at nonzero density. The corresponding relativistic equation 
has been derived in [21] . The pair wave function has the following asymptotic 
form 



where the binding energy is given by e&/2 = m — fi. If we introduce the 
nonrelativistic chemical potential fx — /x — m, we get = —2ft in complete 
agreement with the result obtained by Nozieres and Schmitt-RinkfTo^ (recall 
that ft < in the strong coupling limit). In the nonrelativistic case we may 

D In NJL -type models transition to A ^ phase is accompanied by a sharp increase of 
density towards the value indicated above. Fluctuation corrections should smoother this 
behavior. 



n 



(3) 
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also write // + m = pL + 2m ~ 2m, and then (JIJ) reduces to 

tp ~ ^exp (-r/a a ) , (5) 

where a s = (met,) 1 is the scattering length. This parametrization is very 
convenient for the system of cold fermionic atoms where the scattering length 
is a tunable parameter [21]. However, for quarks the scattering length is an 
ill-defined quantity. 

Let us see, then, what happens to the wave function (j3J) in ultrarelativistic 
case if m — > in line with chiral symmetry restoration. Then 

^-expf-^) 172 )- (6) 



r c V 2 

This result shows that the system develops an instability due to the coales- 
cence of the quark and antiquark branches [21] • The arguments presented 
here are based on the asymptotic form Q of the wave function. To inves- 
tigate the vicinity of the instability point one should consider the system 
of coupled equations for positive and negative energy components written 
down in [21] . This will be done in a forthcoming paper. Therefore the chiral 
symmetry restoration in (NM) — > (QM) transition is a subtle point which 
deserves further investigation. 

Next we turn from the crossover parameter to the Ginzburg-Levanyuk 
number Gi [28j. The two quantities are interrelated since both characterize 
the fluctuating pairs. Let us start with the expression for Gi in case of a 
clean three-dimensional superconductor 



where ((3) ~ 1.2, T c is the critical temperature and Ep - the Fermi energy. 
We note in passing that in [2H] Gi was underestimated by two orders of 
magnitude. This means an extreme narrowness of the fluctuation region. If 
we apply (J7J) to the quark system, take /i for E F and put ji ~ (0.3 — 0.4) 
GeV, T c ~ (0.04 - 0.05) GeV [3El|7j, we obtain 

Gi > 10- 2 (8) 

which is a huge number as compared to that for ordinary superconductors. It 
is by two orders of magnitude larger than the value of Gi for quarks presented 
in |22| . 
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At this point one may ask to what extent is the equation (J7J) for Gi 
applicable to the quark system. One of the ways to derive it [30] is to 
compare the radius of interaction between fluctuations with the correlation 
radius in the system far from the transition point. Then to arrive at (0) use 
is made of the Ginzburg-Landau functional. Therefore the estimate Q relies 
on the applicability of the Ginzburg-Landau theory to the quark system at 
moderate density. We shall discuss this point at the end of the paper. 

Next we wish to express Gi in terms of the crossover parameter. We 
remind the expression for the coherence length 



where vp is the velocity at the Fermi surface. Using Eqs. (J7J) and Q we 
obtain 

We see that the fluctuation effects in a quark system at moderate density 
are very strong. Next we have to determine what is the dominant type of 
fluctuations - that of the order parameter, or of the gauge field potential. 
The fluctuation contribution to the free energy can be estimated j2H] dividing 
the volume of a speciman by the cube of the fluctuation correlation length of 
the given type of fluctuation. The correlation length of the order parameter 
is, as we have seen, ~ 1 fm, while that of the gluon field is ~ 0.2 fm 6 |3Tj . 
Therefore fluctuation of the gauge field dominate. 



3 The Ginzburg-Landau Functional 

The standard framework to consider fluctuations is the Ginzburg-Landau 
functional (GLF). In this section we shall present the derivation of the GLF 
for quarks starting from the general form of the effective action common 
for the NJL, gluon exchange, or instanton models. Considerations presented 
below may be regarded as an extension of our paper [TT] . 
We start with the QCD Euclidean partition function 

Z = J DADijDijexp(-S), (11) 

6 Strictly speaking, within the local Landau- Ginzburg theory we may consider only 
fluctuations with the correlation length larger than the pair size. 
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where 

S = J d^i-i^D, -im + z/i 7 4^ + \ j d A xF; v F; u . (12) 

In (JT2j) color and flavor indices are suppressed, Nf — 2, N c — 3, and the 
chemical potential \x is introduced. Performing integration over the gauge 
fields one gets effective fermion action in terms of cluster expansion ^Hl EH] 

Z = J DipDipexpi- J d 4 xL - S eff ), (13) 

with L = ^(—ij^dfj, — im + ifij^ip and effective action 

S eff = J2n=2^. < ° n >> where = I d i xi){x)g^ tl A c l{x)t a i){x) and the 

double brackets denote irreducible cumulants [32.. 

The derivation of the GLF from the effective action entering into (|T3*j) is 
a complicated task for further work. Only the first step in this direction has 
been done recently ^Hj- Here we shall follow much simpler route which is 
used to derive the GLF for electronic superconductors [BB]- 

We replace S e ff in ()13j) by an effective four- fermion interaction with sym- 
metry properties of the two-flavor QCD. This might be either NJL, or instan- 
ton vertex, or contact interaction imitating one-gluon exchange. Symbolically 
such interaction looks like 

L int = g$kiP)(ijjkiP), (14) 

where the constant g has a dimension m~ 2 . 

Then we perform the Fierz transformation in the Lorentz, color and flavor 
spaces. As a result the interaction term takes the form 

L int = g(ip c RiP)$R + < l p c ). (15) 

Here {ip c Rip) is a scalar diquark in color 3 state 

tp c Rip = iflfiS^e^^ij^pj: (16) 

where C = 7274 and the presence of the color structure 5^6^ signals that 
color symmetry is broken by such diquarks. 

Next step is to integrate the partition function over quark fields. The 
interaction is quartic in the fermion fields and therefore one has to apply the 
Hubbard-Strotonovich trick (bosonization) 

exp(# f dr [ dr(i) c Ri)){^R + ^ c ) 
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DA* DA exp{- dr dr 



\A 



2 
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A{i)R + i) c ) - A*{i) c Ri)) 



(17) 

where A is a complex scalar field. From the Lagrange equation of motion for 
A* we see that A = g(ip c Rip), so that A has a dimension of mass as it should 
be for the gap parameter. Now we can integrate the partition function over 
the quark fields and obtain 

Z = J DA*DAexp^-- dr J dr\A\ 2 - Q B (T, /i, A*, A)| , (18) 

where Qb is the Bogolubov functional 



Q B = - l -tr In ( . aT A T R T ^ IT 74 ) . (19) 



We use the following representation 

o'" 7 )- Hi °-i (20) 

C = 7 274, C~\C = - lf i T , dl = "d^ = -~d,. (21) 

In (119)1 the quark mass is for simplicity set to zero. The important assertion 
we are making is that the gauge field may be introduced at the last step by 
replacing the gradient by the covariant derivative and adding the Yang-Mills 
part of the Lagrangian. We may also recall jBl 13 \H\ that the Fierz trans- 
formation from (fT^jl to (JT5j) results in some extra terms in (|T3|) including 
the chiral condensate (ipaidapdijippj). It has been shown in jTI] that in the 
leading approximation the thermodynamics of the system is determined by 
diquarks. Finally, we wish to note that according to (|18j) the partition func- 
tion is obtained by averaging the Bogolubov functional Qb with the Gaussian 
weight factor |A| 2 /g. 

Next we rearrange (|T9"j) using the identity 

tr\n( B AT ^ + J = tr \n AA T + tr \n{l + A' 1 B(A T )- 1 B + }. (22) 
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Figure 1: Quadratic and quartic terms of GLF. 



Before we apply (|22|) to Qb given by (|19|). we rewrite the operator R (see 
(Jinj)) as R = RoC, where C = 7274 is the charge conjugation operator. Then 
the following formula is easily worked out to be 

A' 1 AR(A T )- 1 A*R + = R R^ ^ A(p)A*(p) ^ G(k)G(p - k), (23) 

p k 

where A = id^j^ — i^ji, and G is the thermal propagator 

G(q) = : = (24) 

q7 + g 4 74 + ^74 q + ^74 

with g 4 = — -n(2n + 1)T, T = Next we return to formula (jTSjl and write 

ZZo-^expf-Cfi-fio)}, (25) 

with (fi — f2 ) being the A ^ part of the thermodynamic potential. Now 
we can finally combine (jl8|) - (jl9)) and (J23|) - (j25j) and somewhat symbolically 
write 

n-n = ^J^dr J dr\A\ 2 - itrlnjl + R R+AA*Gg} . (26) 

Expanding the logarithm in (|26J) in powers of A 2 we arrive at the GLF. 
Following the standard approach we take into account the spatial variation 
of A in the quadratic term and neglect such dependence in the quartic term. 
The diagrams corresponding to the terms proportional to |A| 2 and |A| 4 are 
shown in Fig.l. 
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We start with the quadratic term in (J26j) . It reads 



n» = ! fdr ( W-htrSJtjY.Y. - " AWA ' W ■ (27) 

9 Jo J 2 p fc (A; + 17/74) (p - fc + ^74) 

The trace over flavor and color indices is 

tr f)C R R+ = JV)(iV c - 1) = 4. (28) 

After this trace is taken, formula (|27p still contains ir' over the Lorentz 
indices (7-matrices) . 

The GLF corresponds to the 'soft loop" approximation 

tr'V^ — ~A + 5p 2 . (29) 

£ (A; + ^74) (p - A; + ^74) 

First we compute the static term A. We have 

A = V V / — 1 (30) 

By simple algebra it is easily seen that 

tT (k + i^)(k - ifi l4 ) ~ 2 { kj + (k - fi) 2 + k 2 + (k + /i) 2 } ' (31) 

We assume that one can neglect the contribution of antiparticles. This means 
that only the first term should be kept in the representation ([31)1 . Then we 
perform the Matsubara summation in (|3U|) and obtain 

A = — *th±. (32) 

Here i/p = 2^ 2 /ir 2 is the density of states at the Fermi surface. Due to color 
and flavor degrees of freedom it is four times larger than the corresponding 
BCS theory factor. Parameter A is the cut-off of the four fermion interaction. 
In BCS theory it is the Debye frequency, in the NJL model A ~ 0.8 GeV. As 
we shall see, the value of A is irrelevant as soon as A/2T C ^> 1. Integration 
in (j3~2|l is over £ = k — \i. 
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To proceed further we make an important assumption. In the summation 
over p in (|27|) we integrate over Fourier coefficients with u 7^ 0. In BCS theory 
this leads to small corrections to the GLF coefficients [32] ■ This means that 
we keep in ()27|) only A(p) and A*(p) having u = frequency. Then 

J>(p)A*(p) = p j ^A(p)A*(p) =p J dr|A(r)| 2 . (33) 

Under the same assumption we can integrate over r the first term in ([27)1 . 
Next we expand the integrand in (|32|) in powers of (T c — T)/T c 

\ th Z~L h Z + ^^( ch Z) \ (34) 



£ 2T £ 2T C 2T C 2 V. 2T C 

Upon insertion into ()27j) and with the account of ()28)1 the first term in the 
right-hand side of cancels with the first term in the right-hand side of 
()27|). This is because of the gap equation or equivalently from the definition 
of the critical temperature. Then we perform the integration in ()32|) keeping 
only the last term in ([3*4^ and setting the limits of integration to =Foo since 
A/2T C ~ 10 and the integrand is a rapidly decreasing function. The result 
for A reads 

Now we turn to the coefficient B in (j2HJ). Keeping in mind the assumption 
formulated prior to Eq. ([3~3~j) . we write 7 



dk 1 / d x " 



B = ^tr' > j „ ft , - p- 



2p 2 /3 j^J (2tt) 3 (A; + z/i74 ) V dkj (k-i^ 4 

v f + , 2V- r + °° <x 7 c(3) 



7 2. / (<* T ,2,2 = -^QfiZi^- ( 36 ) 



48/? ' ^7.^ (£ 2 + fc 4 2 ) 2 ' 967^- 

With the coefficients A and -B at hand we can write down the quadratic part 
of the thermodynamic potential 

{l^=(3v F j dr{t|A(r)| 2 + J^|gradA(r)| 2 }, (37) 



7 In fact the above assumption was tacitly used already in Eq. 
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where t = (T — T c ) /T c . 

The derivation of the quartic term essentially repeats that of the quadratic 
term presented above. Therefore we omit the technical details. The result 
reads 

n(4) = .1 ' <ir|A(r)|4 - (38) 

Now we can assemble the pieces together. We remind the two assumptions 
formulated above: (i) non-zero Matsubara are decoupled, and (ii) the gluon 
field enters via the covariant derivative. The result reads 

Q = (3 J drF(T, /i, A*, A, A'), (39) 
|t|A| 2 + ^|A| 4 + 7 |DA| 2 } - \A{V 2 A l k . (40) 



t = v F 

Here 

D = V - in 

2 



D = V - ig^A 1 , (41) 



P = ^L , = ^L (42 ) 

Note that 7 = £ 2 , where £ 2 is the BCS theory correlation length (see Eq. (8)). 
The last term in (|40j) is the contribution of the gluon field in the Coulomb 
gauge with cubic and quartic terms neglected. 

One may argue that expression (|4f)|) contains nothing new. It coincides 
with the standard GLF with the coefficients (}4*2*|) obtained by Gorkov almost 
half a century ago [34J. For the quark system the GLF was presented in 
[TJ ES] and some later papers. What we have tried to do here is to expose 
several subtle points in the derivation of the GLF for quarks. It should be 
clear that even if one takes for granted that in some mysterious way the 
role of the gluon sector reduces to the generation of an effective four-quark 
interaction it is by far not enough to straightforwardly derive the GLF. 



4 Fluctuation Induced Color Diamagnetism 

We have seen that the (NM) — > (QM) transition takes place in the strong 
fluctuation regime. Fluctuating diquarks are the precursors for the Cooper 
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pairs forming the BCS condensate. Below we consider fluctuations of the 
gluon field. This type of fluctuations is responsible for two effects [36j: (i) 
the lowering of the critical temperature (fluctuation diamagnetism), and (ii) 
modification of the order of the phase transition (first order instead of sec- 
ond). Both points with regard to quark system have been discussed in recent 
years [SZl EH] - 

Our emphasis here will be on another aspect of the gluon field fluctu- 
ation phenomena, namely "the emerging phenomenology of ((A a ) 2 )" gluon 
condensate |3U1 |4*U] . In the last few years there has been a growing interest in 
condensates of dimension two |4*T] . It might be that there is some (indirect) 
connection between the phenomenology of (A 2 ) and gluon field fluctuation 
in dense quark matter. There is a further possibility that the problem of 
(A 2 ) is linked to the Anderson transition in quark matter. In order to clearly 
define what we mean by (A 2 ) in the context of dense quark matter we start 
with the derivation of the basic equations. 

The partition function expressed in terms of the variables A, A* and A' 
reads 

Z = J DA* DAD A 1 exp j— /3 J drF(T, /i, A*, A, A 1 ) j = 

DA*DAZ A , (43) 



where F is the GLF (gUJ) and 

DA* exp J drFj. (44) 
Integration over the gauge fields yields 

Z A = exp{-(3u F J dr{t\A\ 2 + ^|A| 4 ) + 

+ ^/(^ ln(2 ^ |A|2 + k2)} - (45) 

Here 7 = VfI- The coupling constant g 2 in (|43jl includes the averaging of the 
kinetic term |DA| 2 over the genetators X 1 . The corresponding calculations 
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may be found in Ref. 02]- Taking the derivative of InZ^ with respect to |A| 
we obtain 

-^-^ Kl \nZ A = u F {t\A\+f3\A\ 3 }+ 
„ f dk 27^ 2 |A| 



The same derivative calculated directly from (jUj) and (|4~UJ) reads 

_ ^^ lnZA = Mt|A|+/3|A|3}+ 
+ 1 g 2 \A\Z A 1 J DA\A 2 )ex V (-P J dvF)} = 

= u F {t\A\ + /3|A| 3 + 7 /| A| (A 2 )}. (47) 
Comparing (|47)J) and (|4*7jl . one obtains 

(A 2 ) = j L>A'(A 2 ) exp{-/3 J drF} = 
m f dk 1 

= 2T j (48) 

where M 2 = 2h>p'yg 2 \ A| 2 . We recognize in M the London penetration depth. 
The average (A 2 ) is the expectation value of A 2 in a fixed bosonic field A. 
Next we integrate (|47|) back over A and obtain 



F = ^{(t + ^ 2 (A 2 ))|A| 2 + ^|A| 4 } 



(49) 



The critical temperature shifted due to gluon field fluctuations corresponds 
to zero of the quadratic term coefficient jHEl EE ESI 

T^ = T C (1-^ 2 (A 2 )). (50) 

We can estimate g 2 £ 2 as g 2 = Areas — 4, £ ~ 2 fm, g 2 C, 2 ~ 400 GeV~ 2 . 
Then the formal absolute upper bound on (A 2 ) following from (J5Uj) would 
be (A 2 ) < (50 MeV) 2 . Later we shall return to the estimate of (A 2 ). 
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The second effect of the gauge field fluctuations is the replacement of 
the second order phase transition by the first order one. This point was 
thoroughly discussed in the literature [HHl E E3 IHH UJ] • A careful look at 
the expression ()49|) for (A 2 ) shows that it induces the |A| 3 term in the GLF. 
To see this consider equation (|4"8j) for (A 2 ) and expand it in terms of M/A, 
where A ~ 0.8 GeV 8 is the cut-off parameter. 

, i2 , TA f vrM M 2 1 TA q , \ i a i 

< A > = ^ t 1 - 2 T + -y - j * * - (51 » 

Then instead of (}4T?|) we have 

F = u F {(t + 9 r^) |A| 2 - f^X 3 |A| 3 + ^|A| 4 } . (52) 

Due to the term proportional to |A| 3 the potential F acquires a second 
minimum at a finite value of |A|. This means that the phase transition is of 
the first order. For further discussion of this question see |37| I38j. 
From (JoTj) we obtain another upper bound on (A 2 ) 

(A 2 ) < ^ ~ (60MeV) 2 (53) 

TV 2 

for T c = 50 MeV, A = 800 MeV. 

We see that the (NM) — > {QM) transition is accompanied by the gluon 
field fluctuations with the typical magnitude (A 2 ) < (50MeV) 2 . What might 
be the meaning of this result? Is this number large or small? The quantity 
(A 2 ) is gauge- variant, but in jlU] it was shown that it attains its minimum 
in the Coulomb gauge. The comparison of the above result with lattice cal- 
culations of (A 2 ) at zero baryon density below and above Tc [SHI shows that 
our value is by more than an order of magnitude smaller than the lattice 
result. One may argue that the two quantities are defined somewhat differ- 
ently (see jSHl^Oj). The physical reason for possible mismatch is that the 
increase of the baryon/quark density leads to the suppression of the gluon 
field [TU1EII- However the suppression by a factor of 10 - 20 seems to contra- 
dict the expected magnitude of the effect [TUj. The possible way to resolve 
this difficulty (if it exists) is to return to the general expression for the GLF 
()40|) and reanalyze the loop diagram which gives rise to the coefficient 7 in 
(HO)) and then to the coefficient £ 2 in front of (A 2 ) in (p)Uj) . 

8 According to 37^ A = Tc (in the high density regime). 
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5 A Marginal Remark on Anderson 
Localization 



Anderson localization (AL) jH] is one of the basic concepts of contempo- 
rary condensed matter physics. Recently it was realized |lo"l 128] that this 
phenomena is important for the description of superconductivity in strongly 
disordered systems. 

The explanation of AL can be given in terms of the two-particle Green's 
function ()29|) (the loop diagram in Fig.l). This diagram embedded into 
the disordered background determines the dynamical momentum dependent 
diffusion coefficient which turns zero at the localization edge. In the theory of 
superconductivity the disorder is created by impurities. According to Ioffe- 
Regel criterion [IE] the phase transition to the AL regime takes place at 



k F l < 1, (54) 

where I is the mean free path. Approaching the localization region the coef- 
ficient 7 in the GLF ()40|) (and hence the coefficient £ 2 in (JoT)|) ) undergoes a 
drastic change 45 . In particular, close to mobility edge (condition (jUIJ)) the 
coefficient 7 = £ 2 in (J4*U)) is substituted by jU] 

-(f) 2 '' 

where Dq = ^v F l is the Drude diffusion coefficient. Formula (|55|) is one of 
several possible expressions for 7. The concrete scenario of the evolution 
of 7 as a function of the disorder strength depends upon the values of the 
critical exponent and parameters (jo^|) - see Refs. [2E1 EH] • Using (fo"5|) we 
may illustrate the idea of possible presence of AL in the quark system in 
the transition region and estimate the renormalization of the coefficient 7 in 
GLF. 

The hypothesis that AL may play some role in quark systems was first for- 
mulated in [37] . The role of impurities was attributed to random components 
of the gluon fields. The behavior of the quark Green's function was analyzed 
and it was shown that localization should not be confused with confinement. 
In the framework of this idea let us estimate the renormalization of the GLF 
coefficient 7 according to ()55j). The characteristic scale of the quark mean 
free path is / ~ 1 fm. This might be, e.g., the average distance between 
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instantons in the instanton vacuum picture of QCD ^H] 9 . For I = 1 fm, 
v F = 1, T c = 50 MeV equation fl55J) yields 7 ~ 1.2 fm 2 , while 7 = £ 2 with 
£ = 2 fm gives 7 ~ 4 fm 2 . Approaching the edge of localization the ex- 
pression ()51|) and the estimate (|53J1 are not valid any more since the inverse 
penetration depth M in (151)1 becomes energy dependent through the energy 
dependence of 7 (expression (|55|) is valid only in the immediate vicinity of 
the mobility edge). With the above estimate for 7 we may return to Eq. (|5U|) 
and obtain a new upper bound on (A 2 ) which is (A 2 ) < (100 MeV) 2 . Hence 
the discrepancy between our estimate for (A 2 ) and that given in Ref.[3*^j may 
be at least partly eliminated by the account of AL. 

6 Summary 

It is clear that further work is needed before the process of quark matter 
formation is finally understood. Three important features of this process 
have been examined in the present paper. They are: (i) crossover from strong 
coupling/low density to weak coupling/high density, (ii) strong fluctuation 
regime, (iii) possible Anderson localization. We have seen that these three 
points are interrelated. 

It can be definitely stated that there is no direct transition from nuclear 
matter phase to quark superconducting .BCS-like phase. In between the two 
phases there is a crossover region with strong fluctuations and possible An- 
derson localization. Despite the extensive investigations our understanding 
of the (NM) — > (QM) transition is in no way being complete. In particular, 
the behavior of the string tension in this region is almost intractable (see, 
however, Ref-fTS]). 

We have presented a detailed derivation of the GLF within the mean-field 
three-dimensional effective theory. We have found that the GL approach is 
on a shaky ground in this region. Small correlation length of the gluon field 
makes the local approximation baseless, fluctuation corrections are large and 
the influence of antiquarks is not negligible. Therefore the investigation of 
the interior of neutron stars using GLF (see, e.g., [49j) can hardly bring 
reliable results. 

We have found an interesting intersection between the QM physics and 
Anderson localization in highly disordered media. The behavior of the quark 

9 Because of quark spin-flip instantons may be similar to magnetic impurities. 
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loop diagram shown in Fig.l embedded into the gauge field background 
should be studied more deeply. 

Finally we wish to mention that the location of the crossover region on 
the QCD phase diagram has been discussed in a recent publication [50] . 
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